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Abstract 

Recent work in the literature has studied a version of non-commutative Schwarzschild black holes 
where the effects of non-commutativity are described by a mass function depending on both the 
radial variable r and a non-commutativity parameter 9. The present paper studies the asymptotic 
behaviour of solutions of the zero-rest-mass scalar wave equation in such a modified Schwarzschild 
space-time in a neighbourhood of spatial infinity. The analysis is eventually reduced to finding 
solutions of an inhomogeneous Euler-Poisson-Darboux equation, where the parameter affects 
explicitly the functional form of the source term. Interestingly, for finite values of 6, there is full 
qualitative agreement with general relativity: the conformal singularity at spacelike infinity reduces 
in a considerable way the differentiability class of scalar fields at future null infinity. In the physical 
space-time, this means that the scalar field has an asymptotic behaviour with a fall-off going on 
rather more slowly than in flat space-time. 
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I. INTRODUCTION 



The present paper, devoted to the scalar wave equation in a non-commutative 
Schwarzschild space-time, is strongly motivated by three branches of modern gravitational 
physics: 

(T) In their investigation of quantum amplitudes in black-hole evaporation |l| , the authors of 
[sl have considered emission of scalar radiation in a black-hole collapse problem, assuming 
non-spherical perturbations of the scalar field on the final surface T^p, and that the intrinsic 
three-metric describes an exactly spherically-symmetric spatial gravitational field. 

(ii) In general relativity, unexpected features of the asymptotic structure are already found 
to occur: massless scalar fields which have a Bondi-type expansion in powers of near 
null past infinity do not have such an expansion near future null infinity; solutions which 
have physically reasonable Cauchy data may fail to have Bondi-type expansions near null 
infinity ^. 

(iii) According to the models studied in j^, [gI, [3], the non-commutativity of spacetime 
can be encoded in the commutator of operators corresponding to spacetime coordinates, i.e. 
(the integer D below being even) 

[xf',x'']=ie^'\ = 1,2,...,D (1.1) 

when the antisymmetric matrix 6^" is taken to have a block- diagonal form 

9^" = dmg(9i,...,eD/2) 

with 

e, = el Vz = i,2,...,D/2, (1.2) 

n 

the parameter 6 having dimension of length squared and being constant. As shown in [6|, 
the constancy of 6 is very important to obtain a consistent treatment of Lorentz invariance 
and unitarity. The authors of solve the Einstein equations with mass density of a static, 
spherically symmetric, smeared particle-like gravitational source as (hereafter, in agreement 
with our earlier work 8|, we use G = c = h = 1 units) 



which therefore plays the role of matter source. Their resulting spherically symmetric metric 
is 

2m{r,e)\ 2 , f, 2m{r,'''^ 



ds' = -\^l '-^jdt'+l^l dr' + r\dQ' + sin' Qdip'), (1.4) 

where, in terms of the lower incomplete gamma function 

^(i^) = ['"Vte-Ut, (1.5) 



Bis 
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, 2M /3 r2\ ^ ^ 

m(r,^) = — 7 -,- . (1.6) 



we define the mass function 



0F ' \ 2' 49^ 

Thus, if one tries to study emission of scalar radiation as in but in the presence 

of a non- vanishing 9 parameter (cf Sj for the pure gravity case), one is naturally led to 
study a scalar wave equation in a spherically symmetric spacetime whose metric is affected 
by 9, which is the goal of the present paper. Section 2 builds conformal infinity for the 
space-time with metric (1.4). Section 3, following jj], turns the scalar wave equation into 
an inhomogeneous Euler-Poisson-Darboux equation. Section 4 solves such an equation and 
shows under which conditions there is full qualitative agreement with general relativity. 
Concluding remarks and open problems are presented in section 5, while the appendices 
describe relevant mathematical details. 

II. CONFORMAL INFINITY 



Inspired by general relativity, we define a new radial coordinate r* in such a way that 

dr 

1 - 



= (2.1) 



r 



This equation is solved by 



r* = r + 2 ^^T^.dr, (2.2) 

r — 2m[r, 9) 



and the metric (1.4) can be re-expressed in the form 



r 

where 6' G [0,7r],(y9 G [0,27r], and we have defined the 'retarded' coordinate 

u = t-r* g]-cx),cx)[, (2.4) 



and the 'advanced' coordinate 



v = t + r* e]-oo,oo[. (2.5) 



The equations (2.2), (2.4) and (2.5) yield 
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{v — u) 



r J r-2m(r,e)"' 



(2.6) 



r-2m{r,e) 

To lowest order, Eq. (2.6) is solved by ^ , and on defining 



(v — u) 
' 2 



(2.7) 



one finds, by iterated approximations, the asymptotic expansion 



1_ 2 ^,fM^^(FW,, p^^, 



r (v — u) (v — u)"^ \ (v — u)^ 

The limit v +(yo with u, 0, fixed defines future null infinity the limit u +oo 
with v,Q,(p fixed defines past nuU infinity X~, while the limit u —oo,v +oo with 
{u + v),Q,if fixed defines spacelike infinity, i.e. the point /°. The figures below show the 
behaviour of the denominator y = 1 — ^"'^^'^^ in Eq. (2.1) for various values of 9. Interestingly, 
the occurrence of 9 does not introduce new singularities with respect to general relativity. 

It is actually simpler to introduce the coordinates u and v separately, which yields the 
conformally rescaled, "unphysical" metrics (here / = r^^, and dE^ is the metric on a unit 
two-sphere) 

= f [-(1 - 2mf)du^ - 2du dr + r^rfS^] = -{p - 2mf)du^ + 2du df + dT.\ (2.9) 

and 

dS^ = -{f - 2mf)dv^ - 2dv df + dE\ (2.10) 

These metrics are manifestly regular and analytic on their respective hypersurfaces f — 0, 
since their determinants are equal to — sin^ © for all /, including f — 0. The physical space- 
time corresponds to / > in (2.9), and we can extend the manifold to include given 
when / = 0. Similarly, in (2.10), the physical space-time corresponds to / > and can be 
extended to include T~ , given when / = 0. Only the boundary X = X"*" U X~ is adjoined to 
the space-time. 
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In common with general relativity, we note here a difficulty that is encountered if we try 
to identify X~ with X+. If we do extend the region of definition of (2.9) to include negative 
values of /, and then make the replacement / —f, we see that the metric has the form 
(2.10) (with u in place of v) but with a mass function —m in place of m. Thus, the extension 
across X involves a reversal of the sign of the mass function, which is incornpatible with Eq. 
(1.6) unless we advocate a discontinuity in the derivative of the curvature joj] across X. It is 
therefore not reasonable to identify X"*" with X~. 

To sum up, we have two disjoint boundary hypersurfaces X+ and X~, each of which is a 
cylinder with topology S"^ x R. It is clear from (2.9) and (2.10) that each of X^ is a null 
hypersurface (the induced metric at / = being degenerate). These null hypersurfaces are 
generated by rays (given by 0, = constant, / = 0) whose tangents are normals to the 
hypersurfaces. These rays may be taken to be the R's of the topological product S"^ x R. 
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FIG. 1: Plot of the denominator y = 1 



2m{r,e) 



in Eq. (2.1) when 9 = 10"^. 



r 



5 



1 



Y -1 
-2 

-3 

20 40 60 80 100 

r 

FIG. 2: Plot of the denominator y = l- ^^^^ in Eq. (2.1) when 6 

1 

. 975 
0. 95 
. 925 

Y 

0.9 
. 875 
.85 

20 40 60 80 100 

r 

FIG. 3: Plot of the denominator y = l- in Eq. (2.1) when 9 
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FIG. 4: Plot of the denominator y = l- in Eq. (2.1) when 
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III. INHOMOGENEOUS EULER POISSON DARBOUX EQUATION 



The coordinates [u, v) defined in (2.4) and (2.5) are not the most convenient for discussing 
the limits which define conformal infinity [4:]. We therefore define (cf. jij) a function W0{x) 
by requiring that wq^x = r^^) should be equal to r* in (2.2), i.e. 

dx 



weix = 



x'^{2xm{x ^,6) — 1) 



which implies 



{x) = -w'q{x) 



(3.1) 



(3.2) 



x'^{l — 2xm) 

Equation (3.1) defines a one-parameter family of monotone decreasing C°° functions taking 
values over the whole real line. The monotone decreasing character of wg is proved by 
imposing that 1 — 2xm > 0. This is indeed satisfied for sufficiently small values of 6, so that 
1 — 2xm ^ 1 — 2xM, which is positive provided x < A C°° inverse function therefore 
exists, which makes it possible to define new coordinates a, b according to (cf 



WgiX = a) = 



dx 



Wg{x = h) 



x'^{2xm — 1) 
dx 



V 



t r* 

2 2 ^ y 



x=b 



U 

'2 



t r* 
"2^2"' 



(3.3) 



(3.4) 



x^(2xm — 1) 

where the integrals (3.3) and (3.4) involve the mass function m = m(r = x~^). On defining 
/ = r^^ as in section 2, one finds from (2.2), (3.1), (3.3) and (3.4) that 



Wg{f{a, b)) = r*{a, b) = Wo{x = a) + Wo{x = b). 



(3.5) 



Moreover, from (3.2)-(3.4), the metric (2.3) in the (m, f , B, ip) coordinates takes the following 
form in the (a, 6, 9, v?) coordinates: 



ds^ = 4(1 - 2mf)gi{a)g2{b)da db + f'^dT. 



(3.6) 



having defined 



Ml (a) = m(a-\ 9), gi{a) = a"2(l - 2aMi{a)y^ = ge{a), (3.7) 

M2(6) = m{b-\ 6), g2{b) = b'^l - 2bM,ib))-' = ge{b). (3.8) 
In the analysis of the scalar wave equation \^'ip = 0, we now rescale the scalar field ip 



according to 



(3.9) 



where f2 is a real positive function such that 



n = o, n^k ^ 0, g'''n^,n^k = o on r 



(3.10) 



The 'unphysical' scalar field ip satisfies the conformally invariant wave equation in 4 space- 
time dimensions, i.e. 

□ +1)^=0, (3.11) 

where Q is the D'Alembert wave operator [4'], and R is the scalar curvature in the 'unphys- 
ical' space-time with line element 



dr = n'ds' = - 2mf)gi{a)g2{b)da db + n'f-'dT.\ 



(3.12) 



\ 



On choosing the conformal factor in the form f2 = (a + 6)/, we therefore obtain the metric 
tensor 

/O G{a,h) 

G{a,h) 

(a + 6)2 

\0 {a + hfsm^Q) 



(3.13) 



having introduced 



( M 2M /3 
m(a, o) = —p=^ 



1 



' V2'4^/2(a,6), 
F{a,h) = f\a,h) - 2m{a,h)f\a,h). 

G{a, h) = 2(a + bfgi{a)g2{b)F{a, h). 



(3.14) 
(3.15) 
(3.16) 



By virtue of spherical symmetry, we look for solutions of Eq. (3.11) as a linear combina- 
tion of factorized terms as 



V^,(a,6,e,^) = ^^F,„(e,y.), 
(a + h) 



(3.17) 



where 1^^(05 V') are the spherical harmonics on S'^. Substitution of the ansatz (3.17) into 
Eq. (3.11) gives 

L[x]=Se{a,h)x, (3.18) 
where L is the Euler-Poisson-Darboux operator j4| 

l{l + 1) 



L = 



dadb {a + by 



(3.19) 
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which depends imphcitly on 6 through a, h defined in (3.3), (3.4), while Sq is the ^-dependent 
source term 



gi{a)92{h)F{a,h) 



(a + hf 



+ -{a + hfg^{a)g2{h)F{aMR{0), 





(3.20) 



having denoted by R{0) the scalar curvature in Eq. (B24). Inspired by [4'], we now write 
the solution of Eq. (3.18) as the sum x° + L^^S, where x° is the general solution of the 
homogeneous Euler-Poisson-Darboux equation L[x] = 0, while is an integral operator 
with kernel given by the Riemann-Green function (see appendix) of L lOj: 



Xeia, b) = x°(a,&) 
having defined 




R{a, b] a, b')S0{a', b')xe{0'' ■, b')da db', 



D{a,b) 



As is described in j4|, x°(a, b) has the general form 



D{a, b) = {a\ b' : < a < a' < b' < b}. 



(3.21) 



(3.22) 



X\a,b) = ia + by+' 



dX A{a) ^fdX B{b) 



da J (a + 6)'+i \db J (a + 6)'+i 



(3.23) 



with A and B arbitrary C'+^ functions. Moreover, the Riemann-Green function of the 
operator L defined in (3.19) is obtained from the Legendre polynomial of degree / according 



to 



111 



having defined ^ 



R{a, b] a, b') = Pi{z{a, b; a , b')), 



f ^ > v^ {a-a'){b-b') + {a + b'){a' + b) 

[a + b)[a' + b') 



(3.24) 



(3.25) 



IV. QUALITATIVE ANALYSIS OF THE / = SOLUTION 

Hereafter we consider for simplicity the case / = 0; the comparison with our Ref. {4] is 
then easier, and all main features are already displayed. Strictly, we consider an asymptotic 
characteristic initial-value problem where data are specified on past null infinity for a G [0, ao] 
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and on the outgoing null hypersurface a = qq = constant. If / = 0, it is clear from (3.23) 
that the characteristic data can be set to 1: x^i^^^b) = 1, while the Riemann-Green function 
in (3.24) reduces to 1 1^: 



Ri=o{a, h- a\ h') = Po{z{a, 6; a', b')) = 1. 



(4.1) 



The inhomogeneous wave equation (3.18) with / = can be solved with the help of a 



contraction mapping, i.e. 



i 



X{a, b) = x°(a, b) + J2 X'^ia, = 1 + ^ x"(a, b), 



n=l 



n=l 



where 



X"(a,&) 




Seia', 6')X"' (a', b')da' db' = 0((a + bf 



(4.2) 



(4.3) 



D{a,t 



and the series in (4.2) is known to be uniformly convergent near spacelike infinity in general 
relativity iJ]. Moreover, in general relativity the partial derivative x,a as a ^ and b is 
fixed remains bounded, as well as the partial derivative x,b as 6 and a is fixed. The 
second derivative x,aa, however, diverges near future null infinity, which implies that the 
presence of a conformal singularity at spacelike infinity affects the behaviour of scalar fields 
on future null infinity, reducing by a considerable amount their differentiability class. Such 
a property is proved by exploiting the integral representation of x,aa, i-e. 

cb ^ 



S^aia, b')xia, b') + S{a, b')x,aia, b') 



db'. 



(4.4) 



By insertion of (4.2) into (4.4), and bearing in mind (4.3), one finds that the possible 
singularities of x,aa are ruled by the integrals ^\ 



lQ{a,b) = / S ^a{0'ib')db' . 



h{a,b) 



(a + b')S^a{ci, b')db', 



(4.5) 



(4.6) 



where Jq remains finite as a — 0, whereas Ji displays a logarithmic singularity as a — 0. 

If / = in Sec. 3 we find, for finite values of 6, the counterpart of (4.5) given by the 
integral 



Jo (a, 6) 



Sg^a{a, b')db', 



(4.7) 
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where, from the asymptotic formulae as a — and b is fixed, we find, for all finite values of 

f^r-.^^ (4.8) 

[a + 0) 

m(a, b) ~ / A/te"* dt ~ M, (4.9) 
Jo 

F(a,6) = (/^ - 2mf){a,b) ~ ^^^^^ + ^W)' (4-10) 

(a + oj^ 

(a + o)'^ 



2 36 



(a + 6)3 (a + 6)4 



(4.12) 



and hence 



2M62 

^0 ~ -7 — —r^ as a 0, (4.13) 

(a + 6)3 

in agreement with the analysis in [4| for general relativity. These approximations should be 
abandoned only if 9 is so large that (cf. (4.9)) 

lim^a^ = constant. (4-14) 

Furthermore, the counterpart of (4.6) is given by the integral 

7i(a,6)= [ {a + b')Se,aia,b')db' -2M\og{a) asa^O, (4.15) 
Jo 

again in full agreement with jj]. Note that a more accurate asymptotic expansion of the 
source term would be 

So ^^{1 - 2bM,{b))-\ (4.16) 

but this does not modify the leading terms as a ^ in (4.13) and (4.15). 
The fi gures below show the behaviour of X,ai X,b sud x,aa- 
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FIG. 5: Plot of the solution x(a,6) in (4.2) and (4.3) with n = 1 at small a, when 9 takes finite 
values and b = 0.3. 
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FIG. 6: Plot of the partial derivative x,a when 9 = 10^^ and b = 0.3. Such a derivative is clearly 
bounded, according to the theoretical expectations The identical behaviour is displayed by x,b 
when 9 = 10"'^ and a = 0.3. 
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o.crs- — Q.i 0.15 0.2 n.?s n. -^ 

FIG. 7: Plot of the partial derivative x,aa when 9 = 10~^ and b = 0.3. The logarithmic singularity 
as a — >^ 0+ is clearly displayed, and it occurs at all finite values of 6. 
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V. CONCLUDING REMARKS 



Ever since Penrose 12|] developed a geometrical picture of an isolated system in general 
relativity as a space-time admitting future and past null infinity (with the associated fall- 
off of the metric along null geodesies going off to infinity), there has been always great 
interest in this coordinate-free way of bringing infinity to a 'finite distance' and discussing 
the asymptotic structure of space-time. Moreover, the conceptual revolution brought about 
by non- commutative geometry [l3 . Q, [l5| has led to an assessment of the very concept of 



space-time manifo 
the work in Refs. 
as follows. 



with 'corrections' to it evaluated, for example, along the lines of 
l6|, |7j. Within this broad framework, the contributions of our paper are 



(i) Construction of conformal infinity for the spherically symmetric space-time which incor- 
porates noncommutative-geometry corrections to Schwarzschild space-time. 

(ii) Evaluation of the source term (3.20) in the inhomogeneous Euler-Poisson-Darboux 
equation which describes the scalar wave equation in the unphysical space-time obtained 
after conformal rescaling of the original metric (1.4). 

(iii) Qualitative analysis of the asymptotic characteristic initial-value problem in the / = 
case, finding again the logarithmic singularity as shown in Eq. (4.15). In the original, 
'physical' space-time with metric (1.4), such a singularity implies that the large-r behaviour 
of the scalar field is described by the asymptotic expansion 

ci C2 , rfilog(r) 



i 



+ 0(r~ 



(5.1) 



and therefore the field falls off at large r rather more slowly than in fiat space-time 

(iv) Numerical support for all results in sections 3 and 4, as shown by figures 5-7 at the end 
of section 4. 

Our results are thus an encouraging progress towards a rigorous theory of wavelike phe- 
nomena in noncommutative geometry, along the lines of the conformal-infinity program of 
Penrose for general relativity 9| . Hopefully, the physical applications to isolated gravitating 
systems in a noncommutative framework, and possibly to black-hole evaporation, will also 
become clear in the years to come. 
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APPENDIX A: THE RIEMANN GREEN FUNCTION 

An hyperbolic equation in two independent variables can be always cast in the canonical 
form [lO] 

/ (9^ d d \ 

= + '^(^' y)^ + K^^ y)Tr + y) x{x, y) = H{x, y). (Al) 



\ dxdy dx dy 

One can then use the Riemann integral representation of the solution 



io|. 



For this purpose, 



denoting by the adjoint of the operator L in (Al), which acts according to 

L'^ix] = X,xy - {ax),x - {bx),y + cx, (A2) 

one has to find the Riemann kernel R{x,y; subject to the following conditions {{C,,'']) 
being the coordinates of a point P such that the characteristics through it intersect a curve 
C at points A and B, AP being a segment with constant y, and BP being a segment with 
constant x): 

(i) As a function of x and y, R satisfies the adjoint equation 

^U[^] = 0, (A3) 

(ii) R^^ = bR on AP, i.e. 

RA^^ 2/; ^) ^7) = K^, v)R{x, y, ^, r])ony = r], (A4) 
and R^y = aR on BP, i.e. 

R,yi^, y, ^, V) = a(^, y)Rix, y, ^, r])onx = ^, (A5) 

(iii) R equals 1 at P, i.e. 

i?(e,r/;e,r/) = l. (A6) 
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It is then possible to express the solution of Eq. (Al) in the form 10 1 



X(P) = -[x{A)R{A) + x{B)R{B)] + 



R 



X,y + aR 



^Ry 



X 



AB 

dy 1 + 



^X,x + (hR - \r,^ ) X 




dx 



R{x,y;^,r])H{x,y)dxdy, (A7) 



where f2 is a domain with boundary. 

Note that Eqs. (A4) and (A5) are ordinary differential equations for the Riemann function 
R{x, y; ^, 7]) along the characteristics parallel to the coordinate axes. By virtue of (A6), their 
integration yields 



R{x,r];^,r]) = exp b{\,r])dX, 



(A8) 



RU,y;^,v) = ^W a{X,^)dX 



(A9) 

which are the values of the Riemann kernel R along the characteristics through P. Equation 
(A7) yields instead the solution of Eq. (Al) for arbitrary initial values given along an 
arbitrary non-characteristic curve C, by means of a solution R of the adjoint equation (A3) 
which depends on x, y and two parameters ^, rj. Unlike Xy the Riemann function solves a 
characteristic initial- value problem jlo| . 



APPENDIX B: SCALAR CURVATURE IN THE UNPHYSICAL SPACE-TIME 

In section 3, the evaluation of the scalar curvature R = R{0) in the unphysical space-time 
with metric (3.13) is as follows. With the notation in Eqs. (3.13)-(3.16), the non-vanishing 
connection coefficients turn out to be (no summation over repeated indices) 

V%^ = T'^^ = -{a + h)G-\ (Bl) 

r;^ = r'^^ = -(« + ^)G'-^in2e, (B2) 

= ^ ^G^a, = G ^G^b, (B3) 

r®e = = r®e = t%, = r^^ = v^, = r^,^ = r;, = (a + br\ (B4) 

r®^ = -sine cos e, (bs) 
r^e^ = r;e = cote. (B6) 



16 



The resulting Riemann tensor is evaluated from the general formula in a coordinate basis 



pA pA pA _j_ pa pA pa p-' 

IJ,vp A'Pi'^ MP CK*^ A"^ CK/f 



^a -pA 



(B7) 



Since we evaluate the scalar curvature R = g'^^R^i,, which in our case is equal to (hereafter, 
no summation over repeated indices a or 6 or or ip) 



eei 



we only need the 12 components of Riemann occurring in 



Rab — R"'nnh + R^nhh + R^nfi,h + R 



aab ~ abb ~ aQb ' aipb ' 



i?e0 — -R" 



0a0 "T -"-060 + 000 + -R ©990) 



p pa I p6 I p0 I TDip 



Among these, the only non- vanishing are 



~G ^G^ab + G "^G^aGfi, 



aab 



-"■0<^0 



= l-2G'-\ 
R%e^={l-2G-')s\n^Q, 



and hence 



R = -2G-^ ( G 



a.-^)+2(a + 6)-(l-| 



With the notation in Eqs. (3.14)-(3.16), we find 



2 _^ fi-Lo 



G,aG,b 



a + b) gi F 



92,b 



{a + b) g2 F J ' 



G 



+ G 



2 -^,ab _ FgF^b 

F F2 



(a + 6) 



Therefore Eq. (B15) yields 

R{e) = {a + b)-^ 



FaFb 



,ab 



For example, in general relativity 

^ dF df dwoif)dwo{a) 



df dwo(f) dwo(a) da 
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g,{a){2f-6Mf)F 



(B8) 

(B9) 

(BIO) 
(Bll) 

(B12) 

(B13) 
(B14) 

(B15) 

(B16) 

(BIT) 
(B18) 

(B19) 
(B20) 



by virtue of 
and hence 

Fab - = gi{a){2f - mf),bF = gi{a)g2{h){2 - l2Mf)F\ (B22) 



12M/ 
(a + by 



Our definition of scalar curvature has opposite sign with respect to the work in but of 
course this does not affect the results. 
In section 3, Eq. (B19) yields instead 

R{e) = (B24) 
[a + by 

where m is the mass function defined in Eqs. (1.6) and (3.14). 
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